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Abstract 

The equations of motion for a fully ionized hydrogenic plasma in 
applied coaxial electric and magnetic fields are analyzed, where the 
term for the Hall effect in the generalized Ohm's law equation picks 
up a factor of 1/2 relative to its usual expression. Magnetization of 
the medium is incorporated through the decomposition of the Hall 
term and the inclusion of the magnetization force, which is found to 
equal or exceed the gradient of the scalar pressure. A limit on the 
kinetic pressure obtains which corresponds to the usual limit of unity 
for a certain selection of parameters. Solutions of these equations 
for the free motion of the charges in the case of an infinite column 
with azimuthal symmetry are compared for various prescribed pres- 
sure profiles, where one finds that the profile near the outer edge plays 
an important role in the feasibility of the equilibrium. 
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1 Introduction 



The analysis of the many-body system commonly called a plasma is the mar- 
riage of electrodynamics and fluid mechanics. In the classical regime, one con- 
siders Maxwell's theory of electromagnetic fields [pQ] and Newton's theory of 
dynamics [[2]] as expressed by the moments of the Vlasov equation [[HI HI IS]]- 
For this article, we consider the application of electromagnetic hydrodynam- 
ics (EM-HD) to a fully ionized medium of a single hydrogenic species in the 
neutral fluid limit, with consideration of the convective terms and inclusion of 
the magnetization force. The distinction between the free and bound current 
is maintained formally by isolation of the fluid and drift velocities. Colli- 
sional disruption of the gyromotion is modelled by a simple correction factor 
applied to the species magnetic moment and gyro-momentum. The pressure 
across an infinite, symmetric column of unit radius is prescribed, and the 
equations of motion are solved for the motion of the charges as represented 
by the free and bound momentum and current densities. The solutions for 
the free momentum are sensitive to the pressure profile near the outer edge, 
and a limit on the kinetic pressure obtains from the nonlinear magnetization 
model. 

Here we consider the fluid description of a neutral, Zi — 1 plasma of 
species s G {e,i} with total particle density n = J2 s n s = 2^cb mass density 
Pm = Y, s n s m s = n (m e + rrii), charge density p e = Y. s n s e s = 0, free 
momentum density p m V/ = J2 s n s m s^s = n (m e V e + mjV;), free current 
density J/ = n s e s V s = n e(Vj — V e ), and pressure p = nT = n s T s = 
J2 s Ps = n Q (T e +Ti) for T = k^T. Note that the plasma density n is 1/2 the 
total particle density n, that energy equipartition gives + wj = 3T s /2 
such that v^" = ^/2T s /m s , and that V e ,i are the species fluid or guiding 
center velocities. Addressing in turn the equations for the electromagnetic 
fields, mass, energy, and momentum, we arrive at a model which predicts 
the equilibrium current and plasma flow from the prescribed pressure and 
applied electric and magnetic fields. 
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2 Theoretical model 



2.1 Maxwell's field equations 

The classical theory of electromagnetism has been well established, and suit- 
ably quantized [[6l [7]], its predictions have been experimentally verified to 
the level of parts per billion [[H [9]]. The energy density of a magnetized 
plasma falls between that of everyday experience and that of high-energy 
collisions, and there is no reason to believe a priori that the theory should 
fail for a fully ionized medium. Formally, one may express the equations 

V ■ E = p e /e , V x B = /i J + fi e dE/dt , (1) 
V • B = , VxE = -dB/dt , (2) 

in Lorentz covariant notation as 

= (<W" - d u (d fl A fl ) = , (3) 

Q^v = QXpp, + Q^puX + QvpXn = q ^ ( 4 ) 

where the field tensor is the four-curl of the electromagnetic potential, F^ u = 
d^A v - VA*, and its dual is = e^ al3 F a(3 /2, where & va $ is the permu- 
tation tensor. In natural units and intrinsic form, one writes d *F = J and 
d F = for F = d A, and Maxwell's theory amounts to the relation between 
the source current 3-form and the connection 1-form through the curvature 2- 
form, d *d A = J. Current is conserved as a consequence of gauge invariance 
and Noether's theorem [[TUl E]], d^J^ = 0, just as energy and momentum 
are conserved from space-time translational invariance. Here we consider the 
neutral fluid limit p e — > such that V • E — > 0, which distinguishes the 
EM-HD model from that of resistive magnetohydrodynamics (MHD) using 
the quasineutral approximation. The EM-HD model is not restricted to the 
neutral fluid limit; however, the fluid equations of motion become much more 
complicated for a plasma with non-vanishing space charge density. 



2.2 Mass and energy equations 

The equations for mass and energy conservation are taken as the continuity 
equation with adiabatic closure, applicable by species, 

dn s m s D p s 

— — + V • {n s m s V s ) = , — = , 5 

ot Dt (n q m,r 
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where D/Dt = d/dt + V s • V is the convective derivative and noting that 
the effect of plasma heating through resistive dissipation is not incorporated. 
The degrees of freedom for the mass and energy density carried by a species 
are given by the species particle density and temperature. Here, we account 
for these equations by assuming their degrees of freedom in the form of the 
density and temperature profiles 

n s (r) = n s (0) + [n s (r x ) - n s (0)] (r/r x ) a " s , (6) 

and similarly for T s , for radius < r < r\ and profile exponent a = 2. 



2.3 Momentum equations 

The d/dt — > equilibrium equations of motion, retaining from the stress 
tensor only the gradient of the scalar pressure, read 



n s m s (V s • V) V s + Vp s = n s e s (E + V, x B) + F sk + F 



Ms 



(7) 



where F Ms is the magnetization force to be discussed later and k ^ s for the 
friction term F ei = —F ie = m e v ei i^je representing interspecies collisions. 
From the definitions of the momentum and current densities, 



J/ 



n 



rrii m e 
e — e 



V e 



(8) 



we may exchange the ion and electron velocities {Vj, V e } for the pair {Vf, J/}, 

J/ 



V/ + 



epr, 



-VTii 



(9) 



The sum of Equations ([7j) gives the equilibrium net force balance equation 

C + + V P+ -F M + = J/xB , (10) 
and their difference the generalized Ohm's law equation for steady currents 
C_+Vp_-F M - = 2(n eE-Fei) + [2nQeV / -no(m i -m e )J//p m ]xB , (11) 
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where p± = Pi±p e and similarly for the magnetization force Fm±, an d the 
convective terms C± = n Q [rrii (Vj • V) Vj ± m e (V e • V) V e ] are given by 



C + 



n (rrii + m e ) 
n {rrii - rn e ) 
+ 2nQm e rrii 



(V f • V) V/ + 
(V/ . V) V/ - 
(V/ ■ V) 



m e mi 




> 




e 2 


\pm 


Pm 


m e rrii 




< 




e 2 




Pm. 



J/ 

&Pm \ &Pm 



V V 



(12) 
(13) 
(14) 



With a rearrangement of factors, the generalized Ohm's law equation may 
be put into the form 

[C_ + Vp- - F M _] /2n e = E-?7J / +V / xB-(m i -m e )J / xB/2ep m , (15) 

using the resistivity rj = m e v ei /nQe 2 for interspecies collision rate v ei . The 
final term in Equation (TT5T) has acquired a factor of 1/2 relative to its usual 
expression [[4, 12J], which in our notation upon m e rrii would equal 
Jf x B/noe. The reason for the difference is that the standard deriva- 
tion, with its application of extraneous mass factors to the equations of 
motion, does not respect the unwritten factor of units, as follows: given 
Equations ([7j) for the net species forces F« and F e in SI units of force- 
density, Nt/m 3 , the standard derivation would take the difference equa- 
tion to be m e Fi — miF e followed by m e m^; however, the amount of 
force present did not change, and reinstating the unit factors reveals that 
F_(Nt/m 3 ) = [m e (kg)Fi](Nt/m 3 /kg/m e ) - [mi(kg)F e ](Nt/m 3 /kg/mi). One 
must be careful to ensure that all one's quantities are in the same (SI) units 
before adding or subtracting them to achieve the physical result. Follow- 
ing the standard derivation [[!]], one would compare the numerator factor 
(m e + rrii) /rrii ~~ * 1 with 1 + 1 = 2. That factor is important, for it predicts 
that the effect of the Hall term on the driven current is only half of what is 
commonly taken, and upon neglect of the fluid velocity V/ — > 0, one may 
substitute the remainder of the net force balance Equation ([TO]) for the Hall 
term in the Ohm's law equation, whereupon neglect of the pressure gradients 
Vpi e — > and any residual magnetization force F M ± — > yields 



n eE - F £ 



c , n {mi-m e ) c 

Pm . 



(16) 



5 



giving Jf = crE for DC conductivity a = 1/rj. Reinstating the accelera- 
tion term oc dJf/dt recovers the AC conductivity [|13j] and with vanish- 
ing damping factor a plasma frequency of u p = e^/n (mi + m e )/mim e eo- 
Equilibrium requires that neither the momentum nor the current diverge, 
V • (pmY f) = V • J/ = 0. We next account for gyromotion by the substi- 
tutions Vf — > Vf + V d = V and Jf — > Jf + 3d = J everywhere except 
the friction term F ei representing collisional disruption of the free motion, 
as the collisional disruption of the gyromotion will be considered within the 
nonlinear magnetization model. 

2.4 Plasma magnetization 

Consider coaxial applied electric and magnetic fields Eo and Ho- The elec- 
tric field will drive a free current Jf that in turn creates a magnetic field 
Uf which in conjunction with the applied field H + ttf = H produces a 
magnetization M resulting from the gyromotion of the constituent particles. 
Our approach to the treatment of magnetization most closely follows that 
found in References [[Ml [15], EH]], except that we will be fully decomposing 
the Hall term J x B in terms of H and M, and our logic follows 

E ^ J/^Hj + Ho^M, (17) 

from which one finds the drift velocity and the diamagnetic current from 
the curls of the net gyro-momentum and magnetic dipole moment densities, 
respectively: p m V d = V x L g and J d = V x M. The gyrovector is defined 
by Q s = — e s B s /m s , where the field felt by a particle of species s is B s //x = 
H + M — j2 s and points along b s = h = H/H for jl s the magnetic moment 
of a single particle. The net gyro-momentum is the sum of each species 
contribution, L g = Y2 S fs n sh, and similarly for the net magnetization, M = 
fs n sfis, where the particle gyro- momentum is the cross product of the 
gyroradius r gs = vj~/u s and the perpendicular momentum, l s = r gs xm s vjr = 
2F s uj s /u^, the magnetic dipole moment is jl s = —{Wj~ / B s )b s = —T S 'B S /B^, 
and the factor f s representing collisional disruption of the gyromotion is 
modelled as 

/, = ^ (1 - e-/-) | ^ , J* , (18) 

J u s K >\^uj s /v s forcj s <z/ s , v ; 

where v s = v ss + v s \. is the net species collision rate, and is normalized so that 
fp = J2 s fsPs- Scattering times t ss / = \ jv ss i are calculated using common 
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formulas [[H1HB]] asr ee = Celml /2 Te /2 /n e e A In A ee , r ei = Ce^m^T^/v^i^e 2 In A e 
Tjj = Ce^ml Tf' 2 /niefh\Au, and z/ ie = n e m e h> e i/nimi, for constant C = 
127r 3//2 . With plasma parameter a s = (n s — l)/n s approximately unity for 
sufficient density, the species field is B s //x = H + M fc + a s M s , and using p = 
p//i gives a final magnetization model for M = — M/i = 2n (f e fl e + fifli)/2 
of 

M= fePe fjPi = fp_ , , 

H -Mi- a M e H-M e - a Mi H — aM K ' 

which has solution M/H = (1 - \Jl — Aafp/H 2 ) /2a for < a < 1 and 
/3-limit (3 = 2p/H 2 < l/2af on the ratio of the kinetic to the magnetic 
pressure, which goes to 1/2 for a, f — > 1 and to 1/ f for a — » 1/2 as a — > 0. 
As a; s is in terms of B s , an iterative approach to the collisionality factor may 
be defined by M n (H - aM n ) = f n p, starting at f = 1 with f n = / n (M n _i), 
and for a dense, magnetically confined plasma we find / remains very close 
to unity. The decomposed Hall term then reads, as Jo = V x Ho is nonzero 
only for r > r±, 

J x B//i = [V x (H + M)] x (H + M) , (20) 
= [(H + M) ■ V] (H + M) — V|H + M| 2 /2 , (21) 
= [(H + M) • V] (H + M) - HVH - MVM + VMH .(22) 

2.5 Magnetization force 

The magnetization of the plasma medium also gives rise to a magnetization 
force felt by the dipoles in an external field [ [T5l [19]]. Here we consider a 
generalization of the macroscopic force densities [ |20l l2"Tj ] given by Lorentz 
and Kelvin, F LK = Hq3 x H + /^M • VH, and by Korteweg and Helmholtz, 
F KH = J x B - H HV/i/2, modelled as 

Ym = /ioVM ■ H = —[iqVMH , (23) 

for M = —{M/H)H, noting that its presence ensures that the Hall term of 
Equation (1201) reduces to the correct form Jj-xH = (H • V)H — HVH in the 
free-current limit Jd J/ [[22]]. We consider this force an important effect 
neglected in the usual analysis of plasma equilibrium despite its experimental 
applications in fusion [[23]], magnetic fluids [[2H 123 125]], biophysics [ [27] |28 | 
[23 EU1 EI]], and materials science [[321[33l[Ml[3g[36l[3?J 
and addressed for the case of a stationary equilibrium elsewhere [[22]]- In 



7 



the evaluation of Equations (fTUj) and (TTTj) . we take F M ± = F Mi ±F Me where 
Fms — — HoVM s H. There remains to incorporate any effects arising from 
the intrinsic spin of the particles, requiring a properly quantum mechanical 
treatment of plasma magnetization. 



3 Numerical evaluation 

3.1 Restricted equations 

With restriction to an infinite d/dz — > plasma column of meter radius 
ri = 1 with azimuthal symmetry d/d6 — > 0, equilibrium requires that Vf r = 
Jf r = 0. The applied coaxial fields are taken as E = E z and H = 
H z. Then, the fluid equations of motion, Equations (TTOT) and (1TT1) . for a 
sufficiently dense a — > 1 medium, reduce to a system of four scalar equations 
in four unknowns, {Vjg, Vf z , Jfg, Jf z }- The free current is found from the 9 
and z components of the Ohm's law equation, which in this case gives us 
Jfg = and Jf z = n e 2 E /m e i' e i. The enclosed free current as a function of 
radius If{r) = j Q 2Trr'Jf z (r')dr' then determines the azimuthal magnetic field 

Hf = Hf6 = (If/2irr)6. Five iterations were found sufficient to converge the 
collisionality factor / in the determination of the magnetization M, and the 
drift current found from its curl is 

3 d = V x M = -V x ^Hj =HxV^-^VxH, (24) 

and similarly for the drift momentum pmYd = V x L 9 , from which in (r, 6, z) 
coordinates 

f n rr d M rr 9 M y M\ 

3 ^{ - h ^h'- h '^h- j 'h) • < 25) 

The net magnetic field is given by B = yU (l — M/H)H. We next consider 
the solution of the remaining radial equations for both the free current model 
Vj, valid when p H 2 such that M -C H and Jd «//, and the magnetized 
model V? with V m = V™ + V™. The free current model has the solution 



V? 




V fo = \/(^ t + ^ J /^/)^' ( 26 ) 



fz 2n efi H f 



n (mj - m e ) f 2 _ dp_ 



fo dr 



(rrij — m P ) H n /f . 



8 



and the magnetized model has the more complicated solution 



8 P+ _ IT. . 7 P. _ T.n\ r _ m « m < 72 



y :> = ^{-^- F ^ +J ^- JA )^-^- y s- < 28) 

n (mi-m e ) f 2 m e mj 2 4m e m i /ep m \ 1 

e^p,^ rrii — m e J or 



in 



fz 2n Q eB, 



'"'■-•"■"j.-%J.) + %Vr-VZ. (30) 



2ep m \ B 8 J B, 
3.2 Comparative solutions 

In the following we consider a plasma with central density no(0) = 9x 10 19 /m 3 
and central electron temperature T e (0) = 3keV immersed in an electric field 
E = 20mV/m and magnetic field B = 2T. We will compare the solutions 
for a central ion temperature 3^ of 3keV and 12keV ("cold" and "hot"), where 
the outer ion temperature is set equal to that of the electrons, fj(ri) = T e (n), 
first for a pedestal ratio \n = ^o( r i)/"-o(0) = 1/10 and then for a ratio of 
1/vlO (and similarly for T e ). The pressure pedestal ratio Xp S oes as t ne 
product of the density and temperature pedestal ratios. For the first so- 
lution with cold ions and Xp — 1/100, Figure [lj we find that the central 
(3q ~ 5.4% remains well below the limit of 50% and that the diamagnetic 
current is on the order of kA/m 2 . The drift current's axial component serves 
both to suppress and enhance the net axial current relative to the free axial 
current, and the total free current for this configuration is approximately 
2.7MA. The net magnetization force is found to equal or exceed the pressure 
gradient force across the profile, and the net axial magnetic field reflects the 
diamagnetic contribution. The net fluid velocity is on the order of Mm/s, 
with the maximum axial component exceeding the maximum azimuthal com- 
ponent by a ration of 4/1, and the drift velocity is on the order of mm/s. 
The axial fluid velocity for the magnetized model is slightly suppressed near 
the core compared to the free current model. 

For the hot ions with the same Xpi Figure [U while the kinetic pressure 
is well within the limit, (3q ~ 13.6%, the solution for the azimuthal fluid 
velocity, both free current and magnetized model, has been driven complex 
near the outer edge of the column, taking the axial velocity along with it— 
this equilibrium is unfeasible. Increasing the outer pressure so that Xp = 1/10 
alleviates the difficulty for both the cold and hot ion configurations, Figures [3] 
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and H] respectively, with a total free current about 3.4MA. Reducing the 
applied magnetic field to B Q = IT for this Xp returns a feasible solution 
for the cold ion configuration with /3q ~ 21.7%, Figure [5], which displays a 
marked reduction in the axial fluid velocity (but not the the azimuthal fluid 
velocity) relative to the higher field case. For this field strength, the hot ions 
are simply too hot, with a central (3q ~ 54.4% in excess of the limit. The 
magnetization has been driven complex in the plasma core, affecting every 
quantity dependent upon it and leading to an unphysical solution. Note 
that the quantities pertinent to the free current model are not affected by 
this calamity, indicating that the models may be distinguished for a suitably 
designed experiment. 

4 Conclusions and outlook 

The equilibrium equations of motion for an axially symmetric, magnetized, 
hydrogenic plasma column in applied electric and magnetic fields in the neu- 
tral fluid limit are investigated. Respecting the units for the net difference 
of the species forces results in a factor of 1/2 on the Hall term in the Ohm's 
law equation relative to its usual value. The equations of motion, including 
the macroscopic magnetization force, are used to determine the free momen- 
tum and current densities from prescribed species density and temperature 
profiles, and Maxwell's theory is used to determine the electromagnetic fields 
from the source charge- current density in the neutral fluid limit. The dia- 
magnetic current and drift velocity are found from the net magnetization and 
gyro-momentum using a nonlinear magnetization model including a collision- 
ality correction factor /. The ratio of kinetic to (free) magnetic pressure is 
found to be limited by a (3\{ m which ranges from 1/ f > 1 for a sparse plasma 
to 1/2/ w 1/2 for a dense, magnetically confined plasma. 

The restricted equations yield analytic solutions for the fluid flow and 
current which are computed for a variety of parameter profiles for both the 
free current and fully magnetized models. The pressure pedestal ratio is 
found to affect the feasibility of the equilibrium through the fluid velocity 
near the outer edge of the column, and comparison of the solutions for cold 
Ti = T e and hot % > T e ion temperature profiles indicates that the attainable 
central pressure is given by the limit on the plasma magnetization. The 
free current model neglecting magnetization is distinguished by having no 
pressure limit. 
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The EM-HD model given above is ripe for extension in several different di- 
rections. A primary difficulty to overcome is the reduction of the stress tensor 
to the gradient of the scalar pressure, rather than a gyrotropic or gyroviscous 
tensor more appropriate for a magnetized medium of free charges [[HI US]]. 
Reinstating the acceleration terms in the equations of motion would yield 
the EM-HD dielectric tensor, and incorporation of a dielectric or conductive 
boundary material would yield a model more descriptive of actual devices of 
fusion [[16]] and propulsion [ [4"?t 08]] interest. Evaluating the theory for a 
V = (d/dr,d/d9,d/d(f> — > 0) geometry corresponds to a model for a toka- 
mak, and reallowing d/d<f) could be used for stellerator analysis. Lifting 
the restrictions of the neutral fluid limit requires readdressing the force bal- 
ance equations for non-vanishing space charge density, with an aim towards 
a manifestly covariant description of Maxwell-Minkowsky electrodynamics 
appropriate for a fluid of free charges. 

5 Appendix 

The resistive magnetohydrodynamic (MHD) equations as usually defined in 
the quasineutral approximation refer to a system of 14 scalar equations in 
14 scalar variables, hence are determined to be complete and soluble. These 
equations are a combination of Navier-Stokes and a subset of Maxwell's. 
However, one of the vector equations is actually an identity when viewed 
from the potential formulation of electrodynamics, hence does not determine 
any degrees of freedom. Only by reinstating Gauss's law does the system of 
equations become closed, allowing for the prediction of both the current and 
momentum from the equations of motion. 

Many authors [H El SSI EDI EH El E31 [SI]] define the low frequency 
resistive MHD equations as the zeroth and first order moments of the Vlasov 
equation with adiabatic closure in conjunction with the two curl equations 
among Maxwell's. For the neutral fluid, the sum and difference of the ion 
and electron equations of motion give the net force balance equation and the 
generalized Ohm's law. Using D/Dt = d/dt + Vj • V, we write the usual 
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equations (which do not distinguish between Jf and 3 a nor H and M): 

^ + V.(p m V /H 0, ^ = 0, (31) 

p m ^ = JxB-Vp, V 3 = E + V f xB , (32) 

V x E = , V x B = /i J , (33) 

where is the resistivity and 7 is the appropriate index for the case un- 
der consideration, and the degrees of freedom are pressure p, mass density 
p m , flow velocity V/, current J, and electromagnetic fields E and B, giving 
a naive counting of 14 scalar equations for 14 scalar variables. However, 
while for decades [[HI [55]] the argument has been made that Gauss's law 
may be neglected with impunity, no one within the plasma physics com- 
munity has denied the applicability of the potential formulation of electro- 
dynamics [pQ]. The unnamed of Maxwell's equations (often called the "no- 
monopole" equation, VB = 0) is brought into play during the determination 
of plasma equilibrium, via solution of the Grad-Shafranov equation [[56[ 157]] 
in toroidal geometry or otherwise, which by the naive counting of above 
would introduce an additional scalar equation, thus over-determining the 
system, yet is commonly known simply to allow for the expression of the 
magnetic field in terms of the vector potential, B = V x A. The reason 
doing so is valid is because vector identities by mathematical definition do 
not determine any degrees of freedom; they reduce them. Inserting that ex- 
pression into Faraday's law [[15]], we recover V x (E + dA/dt) = 0, whence 
E = — V$ — dA/dt, which clearly displays the division of the electric field 
into static and dynamic components and reduces three of our naive degrees 
of freedom down to one for which we have no equation. Unless one wishes 
to invent new physics, the resolution is clear — the reinstatement of Gauss's 
law, V • E = — V 2 $ — <9(V • A)/dt = p e /e which vanishes for a neutral 
fluid, is required to close the system of equations, bringing the number of 
scalar equations and degrees of freedom into agreement with the number 
14 — 3 + 1 = 14 — 2 = 12. We remark that Faraday's law is no less an identity 
than the no-monopole equation as both are given by the general theory of 
vector fields. Gauge invariance plays a special role in the local conservation 
of charge, best expressed in manifestly Lorentz covariant notation. 

From a particle physicist's field-theoretic point of view [[HJ [?[ [58l [59l 60]], 
the Maxwell field tensor F^ u = d^A u — d v in media is known to have only 
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3 physical degrees of freedom embodied by the four-potential = ($/c, A) 
subject to the gauge condition, not 3 for each of the electric and magnetic 
fields, which couple to sources given by the conserved four-current = 
(cp e , J) through the inhomogeneous Maxwell equations d^F^ = (d^d^)A u — 
d u (d fl A fl ) = fioJ u , which are explicitly Lorentz covariant and also gauge in- 
variant, and the homogeneous Maxwell equations, given by the divergence 
of the dual tensor F^ u = e fJ,ua ^F a p/2, where e^" 01 " is the permutation tensor, 
as d^F^ = 0, are satisfied identically when written in terms of the electro- 
magnetic potential, hence do not determine any degrees of freedom. Anti- 
symmetry in F^ u immediately implies conservation of the current, dpd^F^ = 
\i§d v J v = 0, thus it carries only 3 degrees of freedom also. One may recast the 
Maxwell equations into a component-free form through the use of differential 
geometry [[?]], where "the existence of integrals implies a duality between 
forms and chains" which may be exploited. In natural units fi = e = c = 1 
and using the exterior derivative d, the Hodge dual *, the connection 1-form 
A = A^dx 1 *, the curvature 2-form F = (—F tiu /2)dx 11 A dx u , and the current 
3-form J = (J x dy Adz + J y dz Adx + J z dx A dy) A dt — p e dx Ady Adz which 
satisfies the continuity equation d J = 0, one writes the field equation as 
d*F = J and the Bianchi identity, which is a statement on the structure 
of the manifold, as dF = 0, whence F = dA, and we remark that gauge 
invariance, through Noether's theorem [[10[ HI]], implies conservation of the 
covariant current. What all this shows is that the natural, physical divi- 
sion of the Maxwell equations is not into the divergence and curl equations 
but rather into the homogeneous and inhomogeneous equations, whereby the 
Bianchi identity carries the structure for the potential formulation and the 
field equation carries the dynamics obtained from the action. 

The implication for plasma physics is clear: the quasineutral approxima- 
tion does worse than just neglect an effect, as it introduces inconsistency into 
the equations when the components of the electrostatic field are treated in 
isolation [[EE]] • Arguing that Maxwell's divergence equations are initial con- 
ditions for the curl equations is incorrect in media, for while in vacuum such 
statement leads to the propagation of electromagnetic radiation with two 
physical states of polarization, the source terms spoil such interpretation, 
and the divergence of the Maxwell- Ampere equation only recovers the equa- 
tion for local charge conservation, which must be respected independently 
of the conservation of mass addressed by the zeroth moment of the Vlasov 
equation, when Gauss's law retains its intended form. Note that authors in- 
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eluding the no-monopole equation explicitly within the system do not make 
the argument of having 14 equations and degrees of freedom, as that equa- 
tion represents an additional member. Claiming that in general the sources 
may be uniquely determined from expressions for the fields is inappropri- 
ate, for while suitable boundary conditions must be supplied, the differential 
operators hence the boundary conditions are applied to the fields, not the 
sources. The reason for the expression "Maxwell-Lorentz electrodynamics" 
is because Maxwell's theory tells one how the fields react to the sources, and 
the Lorentz force through the equations of motion tells the sources how to 
react to the fields; trying to go the other way around the loop is not well 
defined, as the physics is contained within the action from which both the 
field and source equations of motion may be obtained. 

Let us examine in detail where difficulties are encountered by the neo- 
classical approach, a term we use to encompass all non-classical approaches 
to the fluid description of ionized particles regardless of geometry — such dis- 
cussion [[621 E3]] invariably engenders a hostile response [[Ml E5]] from its 
adherents yet is necessary if one is to consider the application of electrody- 
namic field theory in tensor notation to the many-body system commonly 
called a plasma. The scalar degrees of freedom p m and p may be associated 
with the scalar equations for mass and energy conservation, Equations (|3~T|) . 
as no other quantities appear in those equations for the case of vanishing flow 
velocity; the presence of a flow velocity Vy couples those equations to the rest 
of the system to be solved simultaneously. Note that the previous argument 
tacitly assumed that the equations of motion in the form of the generalized 
Ohm's law and the convective force balance, Equations fl32|) . were associated 
with the degrees of freedom {V/, J}; whereas here, without Gauss's law, one 
must determine the electric field from an equation of motion, usually the 
generalized Ohm's law (however the ion [[HE]] and electron [[IS]] equations of 
motion are also used), giving the solution E neo = 77 J — V/ x B. Faraday's law 
in conjunction with the no-monopole equation then relates the electric field 
to the potentials — E raeo = dA/dt + V$, where without Poisson's equation or 
its gauge invariant generalization the relation between the potentials and the 
space charge density p e remains unspecified (in essence, Faraday's law here 
determines a potential $ which is not an independent degree of freedom), 
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and its divergence gives in various gauges 



V ■ (V/ x B — 77J) 

Coloumb (V ■ A = 0) 

d 

Lorenz (V • A = -fi e —^) 
Weyl ($ = 0) 



|-V- A + V 2 $ , 

at 

V 2 $ , 
□ 2 $ , 



(34) 
(35) 
(36) 

(37) 



where the LHS is explicitly gauge invariant whereas the form and interpreta- 
tion of the RHS is dependent upon one's choice of gauge. The issue of gauge 
invariance is a red herring in the discussion, for while true physics must be 
equally described in any and all gauges, the crucial error in the neoclassical 
approach is its use of an equation of motion to determine the electric field, 
which does not respect Lorentz covariance. (Note that modern power gener- 
ators and electric motors certainly are not moving materially at relativistic 
speeds yet make full and practical use of the covariant transformation proper- 
ties of the field tensor through Faraday's law of induction.) Returning to the 
expression for E neo , let us now examine its transformation properties under 
a change of reference frame. Let S be the frame of the neoclassical observer, 
and let S' be the frame moving with velocity Vj with respect to S. With- 
out loss of generality, the flow velocity in S is taken along the thus 
v / = (V/, 0, 0) ^ gives E neo = {rjJ x ,r]J y + V f B z ,i]J z - V f B y ), using Ein- 



stein's velocity addition rule [[HZ]] gives V'j = 0, and for 7 = 1/y 1 — Vj /c 2 

the transformation for proper velocity applies to J, the spatial part of the 
(0, J), giving 



four- current J£ eo 









r]jx 


Keo = r/J' = 


V'Jy 















E' 



(38) 



where E' = [E x , r y(E y — VfB z ), r y(E z + VfB y )] is the transformation law for 
the physical electric field. Equality could hold only if 7/ = 77/7 = 777 implying 
7 = 1, which holds only when Vf = 0, thus only in the neoclassical frame 
of reference but also implying a vanishing flow velocity, contradicting the 
initial assumption Vf ^ 0. The expression for E neo has inherited the nature 
of a velocity vector from its neoclassical determination hence cannot possibly 
represent a true electric field, which does not transform as the spatial part of a 
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four-vector [PQHS1EZ]]- Furthermore, as ultimately B(J) may be determined 
from Ampere's law or the equivalent Biot-Savart law (for steady currents 
only as all the terms with E need be present for Maxwell's theory to respect 
local charge conservation), the neoclassical electric field depends explicitly 
on the two vectorial quantities of current and mass flow, E neo (J , V /). In 
order to completely determine the system, both of those quantities must find 
solution; however, having already used one of our equations of motion in 
the guise of Ohm's law, we have left only one vector equation for the net 
conservation of momentum, p m DVf/Dt + Vp = J x B, which leaves one 
vector's worth of degrees of freedom without solution, leading to the use 
of a stationary equilibrium equation Vp = J x B in the analysis of non- 
stationary plasma experiments [[151 EE]]- We note that the predictions of the 
neoclassical (NCLASS) model for the poloidal velocity found in a tokamak 
presented in Reference [ |66j ] explicitly fail to agree with the experimental 
measurements. By reinstating the determination of the electrostatic field 
via Gauss's law, what returns is the generalized Ohm's law, an equation of 
motion which one may solve for the motion appearing in that equation, which 
in conjunction with the convective force balance equation fully determines 
the system. Ultimately, the various arguments presented in support of the 
neglect of Gauss's law are superseded by the rigorous formalism of differential 
geometry, whereby casting the Maxwell equations into intrinsic, geometric 
form, d*dA = J, comprises very deep and powerful statements concerning 
what is known about our Universe. 
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Figure 1: (Color online.) Solution profiles for \ v = 1/10, n = 9 x 10 19 /m 3 , 
T e = 3keV, % = 3keV, E = 20mV/m, B = 2T, and (3 = 5.4%. (a) Electron 
pressure + and ion pressure x. (b) Kinetic to magnetic pressure ratio (3 + 
and its limit (3u m x. (c) Diamagnetic current in azimuthal + and axial x 
directions, (d) Axial free current + and net current x. (e) Gradient forces 
for pressure — Vp + and magnetization Fm x. (f) Net magnetic field in the 
azimuthal + and axial x directions, (g) Azimuthal fluid velocity for the free 
current model + and the magnetized model x. (h) Axial fluid velocity for 
the free current model + and the magnetized model x . 
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Figure 2: (Color online.) Solution profiles for Xp — 1/10, n = 9 x 10 19 /m 3 , 
T e = 3keV, Ti = 12keV, E = 20mV/m, B = 2T, and (3 = 13.6%. (a) 
Electron pressure + and ion pressure x. (b) Kinetic to magnetic pressure 
ratio j3 + and its limit (3n m x. (c) Diamagnetic current in azimuthal + and 
axial x directions, (d) Axial free current + and net current x . (e) Gradient 
forces for pressure — Vp + and magnetization Fm x . (f ) Net magnetic field in 
the azimuthal + and axial x directions, (g) Azimuthal fluid velocity for the 
free current model + and the magnetized model x . (h) Axial fluid velocity 
for the free current model + and the magnetized model x. Note that the 
solutions for the fluid velocity have been driven complex near the outer edge 
of the plasma column. 
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Figure 3: (Color online.) Solution profiles for Xp — V x/lO, = 9 x 10 19 /m 3 , 
T e = 3keV, % = 3keV, E = 20mV/m, B = 2T, and /3 = 5.4%. (a) Electron 
pressure + and ion pressure x. (b) Kinetic to magnetic pressure ratio (3 + 
and its limit (3u m x. (c) Diamagnetic current in azimuthal + and axial x 
directions, (d) Axial free current + and net current x. (e) Gradient forces 
for pressure — Vp + and magnetization Fm x. (f) Net magnetic field in the 
azimuthal + and axial x directions, (g) Azimuthal fluid velocity for the free 
current model + and the magnetized model x. (h) Axial fluid velocity for 
the free current model + and the magnetized model x . 



25 




Figure 4: (Color online.) Solution profiles for Xp — 1/ \/l0 5 n = 9 x 10 19 /m 3 , 
T e = 3keV, f - = f2keV, £ = 20mV/m, B = 2T, and (3 = 13.6%. (a) 
Electron pressure + and ion pressure x. (b) Kinetic to magnetic pressure 
ratio (3 + and its limit (3n m x. (c) Diamagnetic current in azimuthal + and 
axial x directions, (d) Axial free current + and net current x . (e) Gradient 
forces for pressure — Vp + and magnetization Fm x . (f ) Net magnetic field in 
the azimuthal + and axial x directions, (g) Azimuthal fluid velocity for the 
free current model + and the magnetized model x . (h) Axial fluid velocity 
for the free current model + and the magnetized model x. Note that the 
slight increase in pressure at the outer edge alleviates the difficulty with the 
fluid velocity profile. 
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Figure 5: (Color online.) Solution profiles for Xp — 1/ \/l0 5 n = 9 x 10 19 /m 3 , 
T e = 3keV, f - = 3keV, £ = 20m V/m, 5 = IT, and (3 = 21.7%. (a) 
Electron pressure + and ion pressure x. (b) Kinetic to magnetic pressure 
ratio (3 + and its limit (3n m x. (c) Diamagnetic current in azimuthal + and 
axial x directions, (d) Axial free current + and net current x . (e) Gradient 
forces for pressure — Vp + and magnetization Fm X . (f ) Net magnetic field in 
the azimuthal + and axial x directions, (g) Azimuthal fluid velocity for the 
free current model + and the magnetized model x . (h) Axial fluid velocity 
for the free current model + and the magnetized model x. Note that the 
axial fluid velocity has decreased by a similar factor of 2 compared to the 
B = 2T profile. 
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Figure 6: (Color online.) Solution profiles for Xp — 1/ x/lO, no = 9x 10 19 /m 3 , 
T e = 3keV, f - = f2keV, E = 20mV/m, 5o = 2T, and (3 = 54.4%. (a) 
Electron pressure + and ion pressure x. (b) Kinetic to magnetic pressure 
ratio (3 + and its limit (3n m x. (c) Diamagnetic current in azimuthal + and 
axial x directions, (d) Axial free current + and net current x . (e) Gradient 
forces for pressure — Vp + and magnetization Fm X . (f ) Net magnetic field in 
the azimuthal + and axial x directions, (g) Azimuthal fluid velocity for the 
free current model + and the magnetized model x . (h) Axial fluid velocity 
for the free current model + and the magnetized model x. Note that the 
magnetization has been driven complex near the plasma core, affecting the 
entire magnetization model while leaving the free current model unaffected. 
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